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Abstract 

We prove a shape theorem for the set of infected individuals in a 
spatial epidemic model with 3 states (susceptible-infected-recovered) 
on Z"^, d > 3, when there is no extinction of the infection. For this, 
we derive percolation estimates (using dynamic renormalization tech- 
niques) for a locally dependent random graph in correspondence with 
the epidemic model. 
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1 Introduction 

Mollison has introduced in [13], [E] a stochastic spatial general epidemic 
model on Z'^, describing the evolution of individuals submitted to infection 
by contact contamination of infected neighbors. More precisely, each site of 
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Z*^ can be healthy, infected, or immune. At time 0, there is an infected indi- 
vidual at the origin, and all other sites are occupied by healthy individuals. 
An infected individual emits germs according to a Poisson process, it stays 
infected for a random time, then it recovers and becomes immune to further 
infection. A germ emitted from x G Z'^ goes to one of the neighbors y 
of X chosen at random. If the individual at y is healthy then it becomes 
infected and begins to emit germs; if this individual is infected or immune, 
nothing happens. The germ emission processes and the durations of infec- 
tions of different individuals are mutually independent. 

Since its introduction, this epidemic model has given rise to many studies, 
and to other models that are variations of this "SIR" (Susceptible-Infected- 
Recovered) structure. A first direction is whether the different states aymp- 
totically survive or not, according to the values of the involved parameters 
(e.g. the infection and recovery rates). A second direction is the obtention 
of a shape theorem for the asymptotic behavior of infected individuals, when 
there is no extinction of the infection (throughout this paper, extinction is 
understood as extinction of the infection). 

Kelly in ^ proved that for = 1, extinction is almost sure for the spa- 
tial general epidemic model. Kuulasmaa in [10] has studied the threshold 
behavior of this model in dimension d > 2. He proved that the process 
has a critical infection rate below which extinction is almost certain, and 
above which there is survival. His work (as well as the following ones on 
this model) is based on the analysis of a directed oriented percolation model, 
that he calls a locally dependent random graph, in correspondence with the 
epidemic model. See also the related paper [TT] . 

Cox & Durrett have derived in [3] a shape theorem for the set of infected 
individuals in the spatial general epidemic model on Z^, when there is no 
extinction, and the contamination rule is nearest neighbor. This result was 
extended to a finite range contamination rule by Zhang in [15]. The proofs 
in [5], [15] are based on the correspondence with the locally dependent ran- 
dom graph, and they refer to [5] (which deals with first passage percolation). 
They rely on the introduction of circuits to delimit and control open paths, 
hence cannot be used above dimension d = 2. 

Chabot proved in [3] the shape theorem for the infected individuals of the 
spatial general epidemic model in dimension (i > 3, with the restriction to de- 
terministic durations of infection: in that case the oriented percolation model 
is comparable to a non-oriented Bernoulli percolation model (as noticed in 
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|10j . the case with constant durations of infection is the only one where the 
edges are independent). He also exploited the papers pj by Antal & Pisz- 
tora, and [7j by Grimmett & Marstrand on dynamic renormalization to deal 
with dimension d > 3. He introduced some random neighborhoods for points 
in Z*^ and with these instead of circuits he was able to extend the proof of [5] . 

In the present work, we prove the shape theorem for infected individuals 
in the spatial general epidemic model in dimension d > 3, when the du- 
rations of infection are random variables. There, the comparison with non 
oriented percolation done in is not longer valid. Our approach requires 
to adapt techniques of [7j, and to derive sub-exponential estimates to play 
the role of the exponential estimates of [JLJ. It is then possible to follow the 
skeleton of [3]. 

In Section [2] we define the spatial general epidemic model, the locally depen- 
dent random graph, and we state our main result. Theorem 2A. In Section [3] 



we derive the necessary percolation estimates on the locally dependent ran- 
dom graph for Theorem 2A_ We prove the latter in Section |4| thanks to an 



analysis of the passage times for the epidemic. 



2 The model: definition and result 

Let d > 3. The epidemic model on Z'^ is represented by a Markov process 
{Vt)t>o of state space Q = {0,i,l}^''. The value r]t{x) G {0,i, 1} is the 
state of individual x at time t: state 1 if the individual is healthy (but not 
immune), state i if it is infected, or state if it is immune. To describe how 
the epidemic propagates, we introduce a locally dependent oriented bond 
percolation model on Z*^. 

For X = {xi, ...,Xd)eZ'^,y= {yi, . . . , y^) G Z^, = X;?=i \xi-yi\ 

denotes the norm of x — y, and we write x ^ y if x,y are neighbors, that 
is \\x — 2/11 1 = 1. Let (Tri:,e{x,y) : x,y E Z'^,x ~ y) be independent random 
variables on a probability space, whose probability is denoted by Pa for a 
parameter A > 0, such that 

1) the Tj.'s are positive with a common distribution satisfying P\{T^ = 
0)<1; 

2) the e{x,yys are exponentially distributed with parameter A. 
We define 

X(x,2/) = P ife(x,y)<T.; 
I- otherwise. 

The oriented bond (x, y) is said open with passage time e{x, y) (abbreviated 
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X-open, or open when the parameter is fixed) if X{x, y) = I and dosed (with 
infinite passage time) if X{x, y) = 0. 

For a given infected individual x, denotes the amount of time x stays 
infected; during this time of infection, x emits germs according to a Poisson 
process of parameter 2dX; when is over, x recovers and its state becomes 
forever. An emitted germ from x reaches one of the 2d neighbors of x 
uniformly. If this neighbor is in state 1, it immediately changes to state i, 
and begins to emit germs according to the same rule; if this neighbor is in 
state or z, nothing happens. 

We denote by C° the set of sites s G Z'^ that will ever become infected 
if, at time 0, the origin o = (0, . . . , 0) is the only infected site while all other 
sites are healthy, that is 

C: = {xeZ'':3t> 0,r]t{x) = t\r]oio) = t,^ z ^ o,r/o(^) = !}• (2) 

It was proven in \^ (1-2)] (see also [121 p. 322], [ini Lemma 3.1]) that C° is 
the set of sites that can be reached from the origin following an open path, 
that is a path of open oriented bonds. 

More generally, for each x G Z"^ we define the ingoing and outgoing clus- 
ters to and from x to be 

Ci = {yeZ':y^ x}, = {y E Z'' : x y}, (3) 

and the corresponding critical values to be 

\l = inf{A : Pi\Cl\ = +oo) > 0}, = inf{A : P(|C°| = +oo) > 0}, 

(4) 

where "x — )■ y" means that there exists (at least) an open path T^^y = 
{xq = x,xi, . . . ,Xn = y) from x to y, and 1^41 denotes the cardinality of a set 
A. Although we are using the symbol o for both the origin and the outgoing 
cluster we believe no confusion will arise because in the former case it appears 
as a subindex while in the latter it does so as a superindex. 
We will prove in Section [3] the following proposition. 

Proposition 2.1 

K = K- 

This common value will be denoted by Ac = Xc{Z'^). 
We can now state our main result: 
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Theorem 2.1 Assume X > Ac. Define, for t > 0, 

= {x El/ : X is immune at time t} = {x E Z,'^ : r]t{x) = 0}; 
C,t = {x : X is infected at time t} = {x G Z'^ : ?7f(x) = i}. 

Then there exists a convex subset D eT/ such that, for all e > we have 

(^{l-e)tDr]C°^ C (^UG) C (^(l + e)tDnC°) a.s. fort large enough; (5) 

and if E{T^) < oo we also have 

Ct C ^(1 + e)tD \ (1 — e)tD^ a.s. for t large enough. (6) 

In other words, the epidemic's progression follows linearly the boundary of a 
convex set. 

To derive this theorem we follow some of the fundamental steps of |5] , but 
since in dimensions three or higher, circuits are not useful as in dimension 2, 
this is not a straightforward adaptation. On the percolation model, we first 
construct for each site z E a. neighborhood V{z) in such a way that two 
neighborhoods are always connected by open paths. For x,y E Z*^, we show 
that the time T{x,y) for the epidemic to go from x to ?/ is 'comparable' (in 
a sense to be precised later on) to the time r(x, y) it takes to go from V(x) 
to V{y). Then we approximate the passage times between different sites by 
a subadditive process, we derive through Kingman's subbaditive theorem a 
radial limit (p{x) (for all x), which is asymptotically the linear growth speed 
of the epidemic in direction x. We establish that the global propagation 
speed is at most linear in z for r(o, z). Finally we prove an asymptotic shape 



theorem for r(o, from which we deduce Theorem 2.1 



3 Percolation estimates 

In this section we collect some results concerning the locally dependent ran- 
dom graph, that is the oriented percolation model given by the random vari- 
ables {X{x,y),x,y E Z'') introduced in Section |2j Note that although these 
r.v.'s are not independent, the random vectors {X{x,x + ei), . . . ,X{x,x + 
ed),X{x, X — si), . . . , X{x, X — Gd) : X E U^} (where (ei, . . . , e^) denotes the 
canonical basis of Z*^) are i.i.d. This small dependence forces us to explain 
why some results known for independent percolation remain valid in this 
context. 
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Remark 3.1 The function X{x,y) is increasing in the independent random 
variables and —e{x,y). It then follows as in |3 Lemma (2.1)] that the 
r.v. {X{x,y) : x,y E Z'^, y ~ x) satisfy the following property: If U and V 
are bounded increasing functions of the random variables {X{x,y) : x,y & 
Z^,y ~ x), then E{UV) > E{U)E{V). 

For n G N \ {0}, let = [—n,nY, let dBn denote the boundary vertex 
points of Bn, and, for x G B{x, n) = x + 5„. For A,Rc Z'^, "A ^ R" 
means that there exists an open path F^. ^ from some x G A to some y & R. 

Theorem 3.1 Suppose A < X°, then there exists /3o > such that for all 
n> 0, 

Pxio ^ dBn) < expi-/3on). 

This is a special case of |[2| Theorem (3.1)], whose proof can be adapted to 
obtain: 

Theorem 3.2 Suppose A < A^, then there exists /3j > such that for all 
n> 0, 

Px{dBa o) < exp(-An). 



It is worth noting that in the context of our paper, by Remark 3.1 [21 The- 
orem (3.1)] can be proved using the BK inequality instead of the Reimer 
inequality (see P, Theorems (2.12), (2.19)]). 



Theorems 3.1, 3. 2| yield Proposition 2.1 



Proof of proposition 2.1. Suppose A < A*. Then by translation invariance 
and Theorem 3.2 we have that for any x G dBn, P\{o x) < exp(— /3jn). 
Adding over all points of dB^ we get P\{o — ?■ dBn) < K'n'^ exp{—Pin) for 
some constant K', which implies that lim„_!.+oo Px{o — ?► dBn) = 0. Therefore 
A < A° and A* < A°. The other inequality is obtained similarly. □ 

From now on, we assume A > \J(L^) and use the following notation: For 

x.y^U^.Ac Z<^, 

(i) {x — )■ y within A\ is the event on which there exists an open path 
^x,y = (a^o = = y) from x to y such that Xj G A for all 

z G {0, . . . , n — 1}. Note that the end point y may not belong to A. 
(a) {x — )■ y outside A} is the event on which there exists an open path 
^x,y = {xq = x,xi, . . . ,Xn = y) from X to y such that none of the Xj's 
(i G {0, . . . , n}) belongs to A. 
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Definition 3.1 For xeZ'^,AcZ'^ let 



C* (A) = {y E A : y ^ X within A} and 
C°{A) = {y E A: X y within A}. 

Note that with this definition Q.(v4) C A and C°{A) C A. 

Next proposition on percolation on slabs follows from the methods of [7] 
or [ni Chapter 7]. 

Proposition 3.1 For any k eN\ {0}, let Sk = Z'^'^ x {0, 1, . . . , A;} denote 
the slab of thickness k. Then for k large enough we havemf^eSk P)y{\Cx{Sk)\ = 
+oo) > and inf,.e5, Pxi\C°{Sk)\ = +oo) > 0. 

More precisely, to adapt [51 Chapter 7] , it is convenient to define the processes 
for different values of A on the same probability space, whose probability 
is denoted by P. This is done as follows: Let ei {x, y) be a collection of 
independent exponential r.v.'s as before, but now with parameter 1. Then 
let ex{x,y) = \~^ei{x,y), and 

X,{x,y) = \^ ifeA(a;,y)<r.; 
I- otherwise. 

The following lemma implies that given K and 6i > 0, there exists ei > 
such that for any A E [0,K] the random field {Xx+Si{u,v) : u,v E Z'^} is 
stochastically above the random field {ma.x{Xx{u,v),Y{u,v)} : u,v E Z'^} 
where the random variables Y{u, v) are i.i.d. Bernoulli with parameter ei and 
are independent of the random variables Xx{u,v). This lemma justifies the 
use of this sprinkling technique. Its proof is elementary and will be omitted. 



Then, with Lemma 3.1 one can adapt the proof of [HI Theorem (7.2)] to get 



Proposition 3.1 



Lemma 3.1 Let K and 6i be strictly positive, then there exists e > such 
that for any A < K,u E Z'^, 

P{Xx^Siiu,v) = 1 Vf ~ mI Xx{x,y) : x,y E Z'^, x ^ y) > e a.s. 

We introduce now some notation: For A G Z'^ we define the exterior vertex 
boundary Ay A as: 

AyA = {x E Z'^ : X ^ A, X ^ y for some y E A}. 

If X ^ y let D{x,y) be the smallest number of bonds required to build 
an open path from x to y. For A C Z'', x E A,y E AyA, '^D{x,y) < 
m within A" means that there is an open path V^^y using less than m bonds 
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from X to y whose sites are all in A except for y. 

The rest of this section provides some upper bounds for the tail of the con- 
ditional distribution of D{x,y) given the event {x — )■ y}. These estimates 
are not optimal and better results can be obtained by adapting the methods 
of [1]. Instead of getting exponential decays in ||x — y\\i (or in n) we get 
exponential decays in ||a; — (or in n^^'^). We have adopted this approach 
because the weaker results suffice for our purposes and are easier to obtain. 

Lemma 3.2 There exist 5 > 0, Ci > and k eN \ {0} such that 
(z)\fn>0, X e Bn+k \Bn, ye {Bn+k \ Bn) U Ay(5„+fc \ Bn) we have : 

P{x ^ y within Bn+k \Bn) > S. 

(a) Let 

An,m = {z : —k + n < zi < n, — oo < Z2 < k}U 
{z : -k + n < zi < m + k^Q < Z2 < k}U 
{z : m < Zi < m + k, — oo < Z2 < k}. 

Wn < m, Vx G Am^ni'^y e Am,n U AvAm,n, we have: 

P{D{x,y) < Ci{\\x - y\\i + \x2\ + I1/2I) within Am,n) > S. 

Again, the proof of this lemma relies on the methods of [HI Chapter 7] and [7] . 
Since it is not an entirely straightforward adaptation we make some remarks 
that we believe will help the reader. 

Remark 3.2 In Chapter 7], renormalised sites, which we will call r-sites, 
are introduced. These are the hypercubes Bk = [—k, k^ and their translates 
B{2kx,k) = [—k,kY + 2kx, x G Z"^. We will denote by x both the point 
in Z'^ and the r-site centered at 2kx since we believe no confusion will arise 
from this. Loosely speaking, these r-sites are called occupied if they are well 
connected with their neighbors. Crucial for this method is the fact that for 
any given p G (0, 1) the set of occupied r-sites dominates a Bernoulli prod- 
uct measure of density p if k is large enough. In 16, Chapter 7] and j^, 
p is taken above the critical value for non- oriented Bernoulli percolation on 
some subset ofU^, but here it is convenient to take it above the critical value 
of 2- dimensional oriented bond percolation. This choice of p and the corre- 
sponding choice of k guarantee the following property: 

(P) There exists 7 > such that for any pair of r-sites u and v the probability 
that there exists a path of occupied r-sites going from u e7/ to v e7/ which 
uses at most 2\\u — v\\i r-sites is at least 7. 
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From this property (P) one deduces that there exist S > and Ci such that 
for any pair x,y E Z"^ x [—k,kY~^ the probability that there is an open path 
from X to y contained in 71? x [— /c, A;]"'^^ that uses at most Ci||x — y\\ bonds 
is at least 6. 

Remark 3.3 We have to add \x2\ + \y2\ because if x E {z : —k + n < zi < 
n, — oo < 2:2 < 0} and y E {z : m < zi < m + k, —00 < Z2 < 0}, to move 
from X to y staying in Am,n we need to reach first the set {z : —k + n < zi < 
m + k, < Z2 < k} (i.e. to increase the second coordinate until it is positive). 



Lemma 3.3 Let k be given by Lemma 3.2 and let x and y be points in Z . 
For n E Tj let Hn = {x E 7/ : Xi = n} and define the events 

Jn = {x H^^-i-ik within B{x,nk),i = 0, ... [n/2\}n 

{Hy^+i+ik -> y within B{y,nk),i = 0, . . . [n/2j}, 
Gn = {x dB{x,kn), dB{y,kn)^y}, 

where, for any a G M, [aj denotes the greatest integer not greater than a. 
Then, there exists (3 > such that 

P{Jn\Gn) > l-exp(-/3n). 



Proof of lemma 3_^. By translation invar iance we may assume that x is the 
origin. When occurs there are open paths from to -B(o, ik)\B{o, {i — l)k) 
and to y from B{y,ik) \ B{y, {i — l)k) for i = 2, 



,n. 



by part (i) oi Lemma 3.2 



Hence we conclude 

□ 



Lemma 3.4 Let k be given by Lemma 3.2 and let Gn be as in Lemma 3.3 
Then, there exist constants C2, C3 and ^2 > such that for all x,y E 
n e'N we have 



I? 



P{D{x,y) > C2\\x - y\\i + CsinkYl Gn) < exp(-a2n) 



Proof of lemma 3^. Again, by translation invariance we may assume that x 
is the origin and without loss of generality, we also assume that yi > and 



?/2 > 0. By Lemma 3.3 it suffices to show that 



P{D{o,y)>C2\\x-y\\, + C,{nkY\ J„ 
decays exponentially in n. 
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Let r = [?T./2j , and for 1 < j < r let 

Aq = {z : -k < zi < 0, -oo < Z2 <y2 + k}U 

{z : -k < zi < yi + k, 1/2 < ^2 < I/2 + k}U 

{z : yi < Zi <yi + k, -00 < ^2 < I/2 + k}, 
Aj = {z : -(j + 1)A; < 2:1 < -jk, -00 < -22 < Z/2 + ^jU 

{2; : -(j + < 2:1 < 2/1 + (j + 

y2+jA;<Z2<l/2 + (j + l)A;}U 

{2; : ?/i + jk < zi<yi + {j + l)k, -00 < Z2 < y2 + (j + l)k}. (8) 

Note that the sets Aq, . . . ,Ar are disjoint. Figure 1 should help the reader 



to visualize them. On the event J„, we can reach from the origin each of 
these sets by means of an open path contained in Bnk and from each of these 
sets we can reach y by means of an open path contained in B{y, nk). Hence, 
on Jn for each i E {0, . . . ,r} there exists a random point Ui G Bnk H Ai and 
there is an open path from the origin to Ui such that all its sites except Ui 
are in Bnk H (flJ^jAp. If there are many possible values of Ui we choose the 
first one in some arbitrary deterministic order. Similarly, there is a point 
Vi G B{y,nk) fl AAi and an open path from Vi to y such that all its sites 
are in B{y,nk) fl (fl^^jA^). Let Ui and Vi be possible values of U and Vi 
respectively. Then define 




A, 



A, 



A 



Figure 1: the event W3 



Fi{ui,Vi) 
Ei{ui,Vi) 



{Ui = Ui,Vi = Vi}, 

{D{ui,Vi) < C\\ui — Vi\\i within Ai} and 
^u„v, {Fi{ui, Vi) n Ei{ui, Vi)) , 
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where the union is over all possible values of Ui and Vi. 
Now we define a subset of Z"^ 

Ri = Bnk U B{y, nk) U . . . U n U . . . U A^^,^ , (9) 

and we denote by cTj the cr-algebra generated by {T^, e(x, y) : x E Ri, x ~ y}. 
Then, noting that lF^(u.,^.)n*~QlvKj is (Tj-measurable write for i = 1, . . . r: 

Ui,Vi 

Ui,Vi 

> sJ^P [Hu.. Vi) n J„ n {n]r}^w^)) = 5P ( j„ n {n]r}^w^)) , 

where the sums are over all possible values of Ui and V^, the first inequality 
follows from the facts that Ei{ui,Vi) is independent of cTj and both J„ and 
Ei{ui, Vi) are increasing events, the second inequality from part (ii) of Lemma 
|3.2 and the last equality from the fact that J„ is contained in the union of 



the Fj(Mj,t>j)'s which are disjoint. Now, proceeding by induction one gets: 

P [Jn n {rfrlw^)) < (1 - 5yp{Jn). 

Since we can choose C2 and C3 in such a way that the event {Z)(o, y) > 
C2IIX — y\\i + C-iink)'^} does not occur if any of the WiS occurs, the lemma 
follows. □ 

Noting that modifying the constant 02 the statement of the above lemma 
holds for C3 = 1/A;'^, we get: 



Lemma 3.5 (i) Let C2 he as in Lemma 3.4 Then, there exists as > such 
that P{D{x,y) > C2\\x — y\\i + n'^jx y) < exp{—asn) ; 
(11) P{x y\ |C°| = +00, \Ci\ = +00) = 1. 

4 The shape theorem 

Let 

C = {x e Z'^ : X 00 and 00 x}. (10) 
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Remark 4.1 As a consequence of Lemma 3.5 part (ii), if two sites x,y of 
Z*^ belong to C , then x ^ y and y ^ x. 

4.1 Neighborhoods in C 

Definition 4.1 For x E TJ^ , let 

R° = {y E Z'^ : X ^ y outside C} (outgoing root from x); 
R-^ = {y ET/ : y ^ X outside C} (incoming root to x). 

In particular x belongs to R° and Rl^ if and only if x ^ C . 

Our next lemma says that the distribution of the radius of R"^ U -R^ decreases 
exponentially. 

Lemma 4.1 There exists ai = (Ji{X, d) > such that, for all n eN, 
P [{Rl U R},) n dBn ^ 0) < exp(-ain). 



Proof of lemma 4-1 For n G N \ {0}, R° fl dB2n 7^ means that there exists 
an open path — dB2n outside C. This implies that there exists x E dBn 
satisfying o — j- x — )■ dB2n outside C. Similarly, R^ fl dB2n 7^ implies that 
there exists x E dBn satisfying 9i?2n x ^ outside C. Then for such a 
point, either the cluster C° or the cluster C* is finite, and has a radius larger 
than or equal to n. We adapt to our case [HI Theorems (8.18), (8.21)] to get 
the existence of ctq = cro(A, rf) > such that: 

r P(C°n95(x,n) ^0, < +00) <exp(-aon); . . 

\P{CindB{x,n) ^0, |C;| < +00) <exp(-aon). ^ ^ 

Hence 

p((i?°u/?j.)n9S2„^0) 

< P(i?°naS2n7^0)+P(i?;n952n7^0) 

< 2 ^ P{\C°\ < +00, X ^ dB{x,n)) 

+2 P(|Q| < +oo,dB{x ,n) x) 

< 4|9P„| exp(— (Ton) 

which induces the result. □ 
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To define the neighborhood V(x) on C of a site x, we introduce the smallest 
box whose interior contains R° and R]., which contains elements of C, and 
is such that two elements of C in this box are connected by an open path 
which does not exit from a little larger box. For this last condition, which 
will enable to bound the passage time through V(x), we use the parameter 



C2 obtained in Lemma 3.4 



Definition 4.2 Let C = C2d + 2. Let k,{x) be the smallest / G N \ {0} such 
that 

« 9i?(x,/)nJi?°.Ui?^) = 0; 
(m) B{x,l)r]C y^d); 

{Hi) V {y, z) e {B{x, I) n C)^ y z within B{x, C'l). 
Remark 4.2 By (i) above, W^URic B{x, k(x)). 
The random variable k,{x) has a sub-exponential tail: 

Lemma 4.2 There exists a constant a = a{X, d) > such that, for any 
n eN, 

P{k{x) >n)< exp(-(Tn^/'^). 



Proof of lem ma\4.2 We show that the probability that any of the 3 conditions 
in Definition 4.2 is not achieved for n decreases exponentially in n^/'^: 



(i) By translation invariance, we have by Lemma 4.1 



P {dB{x,n) n U ^ 0) < exp(-crin). 
(a) There exist m E N,a2 = cr2(A, d) > such that 

P{B{x,n)r]C = 0) < exp{-a2[n/{m + l)\). 



(12) 



(13) 



Indeed, since {|C*(S'm)| 



-00} and {\C°{S„ 



-00} are increasing 



events, it follows from Proposition 3.1 and the FKG inequality (see Re- 
mark 3.1 ) that for m = m(A, d) large enough we have inf^^Sm Pi\^xi^m)\ = 
\C°{Sm)\ = +00) > 0. Then, (13) follows from the facts that B{x,n) inter- 
sects [n/{m + 1)J disjoint translates of Sm, and events in disjoint slabs are 
independent. 

(Hi) There exists = 0-3 (A, d) > such that 



(14) 



P{3{y,z) e {B{x, n) n Cf, y^z within (5(x, C'n) 
< exp(— (Tsn^/'*). 
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Indeed, if no open path from y to z (both in B{x,n) fl C) is contained in 
B{x, Cn), then D{y, z) > 2{C' — l)n. Given our choice of C this imphes that 



D{y^z) > C2\\y — z\\i + n. Therefore, (14) follows from part (i) of Lemma 



1331 □ 

We define the (site) neighborhood in C of x by 

V{x) = B{x,K{x))r]C. (15) 



Remark 4.3 By condition (ii) in Definition 4-2, V{x) ^ 0. 



By condition (Hi) in Definition 4.2, for all y,z in V(x), there exists at least 
one open path from y to z, denoted by F*^^ contained in -B(x, C k{x)). If there 
are several such paths we choose the first one according to some deterministic 
order. We finally define an "edge" neighborhood r(a;) of x: 



r(a;) = {(y', z) C B{x, k(x)), (?/', z) openjU 
{{y\z')eTl,,y,zeV{x)}. 



Those neighborhoods satisfy 

V{x) C 5(x, k(x)); 



r(a;) C B{x,C'k{x)). 



(16) 



(17) 



4.2 Radial limits 

We now come back to the spatial epidemic model. Indeed, the underlying 
percolation model does not give any information on the time needed by the 
epidemic to cover C. We first define an approximation for the passage time 
of the epidemic, then we prove the existence of radial limits for this approx- 
imation and for the epidemic. We will follow for this the construction in [5]. 

For x,y E Z'^, if x y,x ^ y and Tx,y = {xq = x,xi, . . . ,Xn = y) de- 
notes an open path from x to y, we define the passage time on V^^y to be 
(see @) 

n-l 

r{Tx,y) = '^e{xi,Xi+i) (18) 

i=0 

and, if X = y, we put t(Tx,x) = 0. 
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For x,y E Z'^, we define the passage time from x to y to be 

{inf T{T^^y) a X ^ y,x y, 
0^"'^ ifx = y, (19) 

+00 otherwise. 

where the infimum is over all possible open paths from x to y. By analogy 
with [1], [5], we also define 

T{x,y) = inf T{x',y'); 

x'eV(x),y'£V(y) 

J2 r{y\z') ifr(x)^ 

,z')(iT{x) 

otherwise. 
By Remarks 4.1[ 4.3, T{x,y) is well defined. 



Remark 4.4 // V(x) fl V(y) 7^ 0, i/ien f(x, = 0. 

We now show that if y E C°\ R°, T{x,y) approximates T{x,y). 

Lemma 4.3 For x G Z'^, if y E C° \ R°, we have 

r{x, y) < r(x, y) < u{x) + f(x, y) + u{y). 



Proof of lemma 4-3. Let Tx,y be an open path from x to y such that t{x, y) = 
T{Tx y). Since y ^ R° this path must intersect C. Let Ci and C2 be the first 
and last points we encounter in C when moving from x to y along T^^y By 



condition (i) of Definition 4.2, ci G V(x) and C2 G V(y): indeed (for instance 
for ci), either x E C and Ci = a;, or the point a E dB{x, k{x)) fl F^. ,^ does 
not belong to R° and Ci is the first point on F^. ^ between x and a; we might 
have ci = C2, if V(x) fl V(?/) 7^ 0. We have 

r = F V F V F 

x,y ^ x,ci V J- ci,C2 ^ C2,y 

where V denotes the concatenation of paths, Tx,ci is an open path from x to 
ci contained in B{x, n{x)), Fci,c2 is an open path from ci to C2 and Fc2,y is 
an open path from C2 to y contained in B{y, K,{y))- We then obtain the first 
inequality of the lemma since: 

r{x,y) < t{Tc,,c2) < r{T^^y) = T{x,y). 

To prove the second inequality of the lemma let F^^ be an open path from 
di E V{x) to d2 E V{y) such that r(Fd^^rf2) = T{x,y). Since the open paths 
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Tx,ci from X to Ci and F*^ (which exists by Remark 4.1) from ci to di have 



edges in T{x) (see (16)) the open path T^^^i = Tx.ci V T*^ from x to (ii 
satisfies r(r^^dj < m(x). Similarly, there is an open path Td2,y from d2 to y 
such that T{Td2,y) < ii(y)- We conclude with 



t{x, y) < r(r^,rfj + T{Td^^d2) + r{Td2,y) < u{x) + f (x, y) + m(?/). 



□ 



Lemma 4.4 For a// x,y,z G Z"^, we /iai'e i/ie subadditivity property 

t{x,z) <r{x,y) +u{y) + T{y,z). (21) 



Proof of lemma 4-4 Let Fa^f, be an open path from a G V{x) to b E V{y) 



such that r(x, y) = r(ra^fe). Similarly, let Tc,d be an open path from c G V{y) 
to d G V(z) such that T{y, z) = f{Tc,d) (we might have a = b, c = d or b = c) . 
Since both b and c are in V{y) there exists an open path F^^ from 6 to c such 
that r(r^^) < m(?/) (see Remark 4.1 and (16)). The lemma then follows since 



the concatenation of these three paths is an open path from a point of V(x) 
to a point of V{z) and 



f(x, z) < T(Ta,b) + T{Tb^c) + r{Tc,d) < t{x, y) + u{y) + f(y, z). 



□ 



Before stating our next lemma we introduce some notation. For x, ?/ G Z° 
let _ 

D{x,y) = inf D{x',y'). 

x'£V{x),y'eV{y) 

Note that unlike D{x,y), D{x,y) is always finite. 

Lemma 4.5 There exist constants C4 anc? ^4 > such that 

P(D{x,y) > C4x-y\\i+n) < expi-a^n^^'^), Vx,y G Z'^,n G N. 



Proof of lemma \4^ Let i^T be a positive constant. Then 

P(D{x, y) > K\\x -y\\i + {2d + l)Kn) 

< P{k{x) > n) + P{K{y) > n) 

+ P(D{x, y) > K\\x -y\\i + {2d + l)Kn, k{x) < n, K{y) < n) 

< P{k{x) > n) + P{K{y) > n) 



16 



+ P{D{x',y')>K\\x-y\\i + {2d+l)Kn,x'^y') 

x'(^B{x,n),y'(^B(y,n) 

< P{k{x) > n) + P{K,{y) > n) 

+ P{D{x',y)>K\\x-y\\i + Kn,x-^y'). 

x'(iB{x,n),y'(iB{y,n) 



Taking K = C2 (given in Lemma 3.5) the result follows from Lemmas 3.5 



andSil □ 

Of course, the random variables u{x) and T{x,y) are almost surely finite. 
But we will need a better control of their size, provided by our next lemma. 

Lemma 4.6 For all x,y E JJ^ , r G N \ {0}, u{x) and T{x,y) have a finite 
r-th moment. 



Proof of lemma 4-6 By Lemma 4.2, u{x) is bounded above by a sum of 



passage times e(?/, z) with y and z in the box B{x, Y), where F is a random 
variable whose moments are all finite. By Lemmas |4.2| and |4.5| the same 
happens to T{x,y). Therefore it suffices to show that if G N) is a 

sequence of i.i.d. random variables and is a random variable taking values 
in N, then the moments of YliLi -^i finite if it is the case for both the 

Xj's and N. To prove this write: 



N 



E{\J2x^n = 5^^(ixi + ...+x„rv=„}) 

n=l 

00 

< Y.i^(\^^+---+^r.np{N=n)Y/' 

n=l 
00 

< 5^[E(|Xi| + ... + |X„|)2'-p(iV = n)]i/2 



1=1 n=l 
00 



n=l 

00 



n=l 



where the second line follows from Cauchy-Schwartz' inequality, the factor 
n^^ counts the number of terms in the development of (|Xi| + . . . + |X„|)^'' 
and the constant depends on the distribution of the Xj's. As has all 
its moments finite P{N = n) decreases faster than n~^''~^ and the sum is 
finite. □ 

We now construct a process {'&.) which is subadditive in every direction. 
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and has a.s., by Kingman's Theorem, a radial hmit denoted by We will 
then check that r(o, ■) also has, in every direction, the same radial limit, and 
we will extend this conclusion to t(o, ■) on the set C° of sites that have ever 
been infected. Hence we first prove 

Theorem 4.1 For all z G Z*^, there exists fi{z) G such that almost surely 



hm ^-^=^^iz), 



lim 



n—^+co 



r(o, nz) 



n 



n 



-fi{z) 



-{nzecs} 



0. 



(22) 



(23) 



Proof of theorem ^.i, (i) For all z G Z"^, (m, n) G N^, let 

dz{m, n) = T{mz, nz) + u{nz). 



(24) 



The process ('(^^(m, n)) (^rn,n)eN^ satisfies the hypotheses of Kingman' subaddi- 
tive ergodic theorem (see [121 Theorem VI.2.6]) by (21). Hence there exists 
/i(z) G such that 



lim —{}z{0,n) = fi{z) a.s. and in L^. 



(25) 



Since the random variables {u{z) : z G I/) are identically distributed, it 



follows from Lemma 4.6 and Chebychev's inequality that X^^o -^(^('^^) ^ 
ne) < +00 for all e > 0, so that by Borel-Cantelli's Lemma 



lim 



u{nz) 



n-5-+oo n 



0, a.s. 



(26) 



Thus by (25), (26) we have (22) for all z G Z^ 



(a) Since R° is a.s. finite, if nz G C°, then nz E C° \ R° for n large 



enough. Hence, from Lemma 4.3, for n large enough we have 



TO, nz) 



n 



-/i(z) 



l{n^ecg\i?°} < 



n(o) + u{nz) 



n 



r(o, n^;) 



^i{z) 



and we conclude by ( 26 ) and ( 22 ) . 



□ 
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4.3 Extending /i 



We have proved the existence of a hnear propagation speed in every direction 
of Z"'. However, to derive an asymptotic shape resuh, in particular for the 
approximating passage times (T{x,y),x,y G Z"^), we need to control this 
propagation speed uniformy in all directions. For this we study fi{z),z G Z*^, 
and we follow [1] to construct a Lipschitz, convex and homogeneous function 
ip which extends n to M'^. The asymptotic shape of the epidemic will be given 



by the convex set D defined later on in (50). 
Lemma 4.7 The function g defined on by 

yzeZ", g{z) = E{^,{0,1)) (27) 
has a bary centric extension to M'^. 



Proof of lemma ^.7, By (25), we have a.s., for all z G Z 



7d 



lim '^=ME (^1 = inf »M = (28) 



To do a barycentric extension of g, we decompose [0, l]*^ in simplexes: each 
of them having a unique barycentric decomposition, g will be defined on its 
elements as the barycenter of its values on extremal points. This (arbitrary) 
construction will be translation invariant. 

1. Let M denote the center of [0, l]"^. We define g{M) to be the mean of 
the values of g on the 2'^ elements of Z'^ fl [0, l]*^. 

2. For the center c of each face F of [0, l]'^ (which is a cube of dimension 
c? — 1) we define g{c) to be the mean of the values of g on the 2'^~^ 
elements of Z*^ fl F. We proceed similarly on each sub-cube, up to 
sub-faces of dimension 2. 

3. Now, on sub-faces F of dimension 2, we link the center to the 4 elements 
of Z'^nF, to obtain 4 triangles, or simplexes, of dimension 2. On each of 
them, we define g{x) for each point x to be the barycentric combination 
of the values of g on the 3 extremal points. 

4. We deal with dimension 3 sub-faces by taking barycentric combinations 
between the dimension 2 simplexes and the center of the dimension 3 
sub-face. This way we have decomposed each dimension 3 sub-face into 
4x6 simplexes, on which for each point x we define g{x) in a barycentric 
way. We go on in the same way until the dimension d cube, that is 
[0,1]'^. 
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The function g is continuous on [0, 1]'^. Then, for z = (zi, . . . , Zd) G U^, to 
each d-cube riiLil-^*' + 1] associate the 2'^~^d\ simplexes translated from 
those described in [0, l]'^, and we define as previously g{x) for each point 
X G rif=i[-^«5^* + 1] be the barycentric combination of the values of g on 
d + 1 extremal points. □ 



Following [H Lemma 3.2], we define a sequence of functions {gn)n>o by 

g{nx) 



n 



(29) 



Lemma 4.8 The elements of the sequence {gn)n>o ore Lipschitz functions 
with a common Lipschitz constant denoted by 7*, hence the sequence is 
equicontinuous. 



Proof of lemma \4~8\ It is enough to prove that g is a Lipschitz function. For 
all x,y G Z'^, by subadditivity and symmetry of -i?. on Z'' we get 



9{x) + g{y) > g{x + y), 
g{x + y)+g{-y) > g{x), 

9{-y) = 9{y) < l|y||i^(ei). 



(30) 
(31) 
(32) 



Indeed, 



gix + y) = E{^^+y{0,l)) = E(T{o,n{x^+y))+uinix + y))) 

< E (r(o, nx) + u{nx)) + E {jinx, n{x + y))) + E {u{n{x + y))) 

= gix) + E{T{o,ny)) + E{uio)) 

= 9{x)+g{y) 



where we have used successively (27), (24), (21) and the translation invari- 



ance of the distributions of r and u. Then, writing (30) for x = {x + y) + {~y) 



gives (31). For (32), we first use the symmetry of i). on Z to get g{—y) 



g{y), that we then combine with (30) to write 



9{y) < ^9{\yih) < \yi\9{ei) = \\y\\ig{ei). 



1=1 



1=1 



Therefore by (30), (31), (32) 



\g{x + y)-g{y)\ < g{y)l{g(^^+y)>g(y)}+g{-y)l{g^^+y)<g(^y)} < ||y||i^(ei). (33) 
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If we now take x,y E M*^, from the previous barycentric construction, let 
(a;o = X, . . . ,Xn = X + y) he the sequence of points on the simplexes crossed 
by [x, X + y]. Then 



n— 1 



\g{x + y)- g{x)\ < ^ \9{xk+i) - 9{xk)\. 

k=0 

Thus, since \\y\\i = Ylt=o W^k+i — Xk\\i, we have to show that on a given 
simplex the Lipschitz constant of g does not depend on the simplex. Assum- 
ing now that x, y belong to the same simplex, they are written uniquely as 
a barycentric combination of the d + 1 extremal points {zq, . . . , Zd) of that 
simplex, zq being the center of the cube translated from [0, l]'^ containing the 
simplex. Similarly, each < i < d, is the barycenter of 2*^ extremal points 
(cj, 1 < 2 < 2"^), with coefficients given by the barycentric construction: 



X 



d d 

^aiZi] y -- 

i=0 

f^icr, Zi 
1=1 1=1 



i=0 i=0 i=0 

2"^ 2'^ 



(34) 



1=1 



1=1 



As {zi — zq,1 < i < d) is a basis of the vector space M'^, denoting by ||.||^ the 
i^-norm w.r.t. this basis, there exists a constant 70 > such that 



V2 G 



7o 



kill < \\z\\*i < 7o|kl|i- 



(35) 



Since [0,1]'^ is decomposed in a finite number 2'^~^d\ of simplexes, (35) is 
valid for all these simplexes, for a constant 7 > which is the infimum of all 
the 7o's. We have, using (33), (34), (35), 

d 



^(oi - I3i)g{zi) 

i=0 

d 

<j2\(^i-Pi\ 



i=0 



\9{x)-g{y)\ = y^^{ai - (3i)g{zi) = ^{at - (3,){g{zi) - g{zQ)) 

i=0 

2d 

^{li - Ki){g{ci) - g{co)) 
1=1 

2-^ 

< ^ |ai - /3i| ^ |ii - Ki|||Q||i5f(ei) 

i=0 1=1 

d 

< ^ |ai - (3i\2'^ X 2 X 2d X ^(ei) = 2'^+^dg{ei)\\x - y\\l 



i=0 
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<2^+V^(ei)||x-i/||i 



hence there exists 7* > such that 

Vx,2/gM^ \9{x)-g{y)\<i*\\x-y\\^. (36) 

□ 

Lemma 4.9 The sequence {gn)n>o converges uniformly on each compact sub- 
set of M!^ to a function (p which extends /i to W^, is Lipschitz with Lips- 
chitz constant 7*, convex and homogeneous (that is which satisfies ip{aix) = 
ai(p{x) for all x eM.'^ and ai > 0). 



Proof of lemma 4-9. (i) For x G Z^, 

g{nx) E{^n.{QA)) ^(^x(0,n)) 



gn{x) 



n n n 



Hence by (25) we get 

hm gm{x) = fi{x) V x G Z'^. (37) 

Let now x G Q'^, and 

= mm{k >l,keN:kxe Z'^}. (38) 

Then, gnN^{x) = gijiN^x) / {uNx) converges to fi{Nxx)/Nx as n goes to in- 
finity. To prove the convergence of gm{x) over the whole sequence, write 
m = n{m)N.j. + j{m) where j(m) G {0, . . . , — 1}, so that 

gijnx) g{n{rn)NxX + j{m)x) 

gm[X) = = r-rz ■ 77 7 

m n[m) IS X + J \jn) 

g{n{m)Nxx) n{m)Nx 

X 



n{m)Nx n{m)Nx+j{m) 
g{n{m)NxX + j{m)x) — g{n{m)Nxx) 



n{m)Nx + j{m) 



By ( 36 ) , the second term of the last right hand side above converges to as 
m goes to infinity. Therefore, 

/ X ,. q(n(m)Nxx) n(m)Nx uiNxx) , , 

lim gm{x) = lim , ; x , , \ , = ^ % Vx G Q"'. 
m-5>oo r?i-s>oo n{m)Nx n{m)Nx+j{m) 

(39) 
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It follows from Lemma 4^ and Arzela-Ascoli's Theorem that any subsequence 
of {gm{x))m>o has a further subsequence that converges uniformly on compact 
subsets of M"' to a Lipschitz function with the same Lipschitz constant 7* 
as the gm^s (cf. (36)). Since, by (39), ip{x) must be equal to fi^N^x) / for 
all a: G Q*^ and ip is Lipschitz, the limiting function does not depend on the 
subsequence and the whole sequence {gm)m>o converges uniformy on compact 
subsets of M'^ to ip which extends /i by (37). 



Vx G Q"', lim gn{x) = ip{x). 



(40) 



This implies convergence on M'^, since every subsequence of {gn)n>o has a 
subsequence which converges uniformly on each compact subset of M*^ to a 
continuous function, equal to ip on Q''. 

(a) To prove that (f is homogeneous we start noting that for 2; G Z*^ and 
/c G N we have: 



lim 

n— >+oo 



lim M^ = /fM = ^M (41) 
n n-5-+oo nk k k 



Now let X eQ'^ and recall that (fi{x) = ix{N^x)/N^ = ip{N^x)/N^. Then if n 
is a multiple of N^, we let k = n/N^ G N and write by (41 ), 

ip{Nxx) ip{kNxx) ip{nx) 



ip{x) 



kN^ 



n 



(42) 



V a; G G N. 



Since A^'^, is a multiple of N^x-, (42) implies: 

^{kx) = = = ^-^^^ = k^{x), 

Hence, if r = n/m and x G Q'^ we have: 

(p{rx) = n(p{{l/m)x) = {n/m)(f{x), 
so that !f is homogeneous on Q'^. 

(in) To prove that ip is convex on Q'', take x,y & Q'^ and a G Q fl (0, 1). 
Then let ki, k2 be elements in N such that kia G N, k2X G Z'^ and A;2l/ G Z*^. 
Using subadditivity of g and homogeneity of write: 



ip{ax + (1 — a)y) 



lim 

n— >oo 



g{nax + n(l — a;)y) 
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^. g{nkiak2X + nki{l — a)k2y) 

n-^oo nkik2 

g{nkiak2x) + g{nki{l - a)k2y) 

< lim — 

n-5>oo nkik2 

_ ip{kik2ax) + (p{kik2{l - a)y) 
kik2 

= aip{x) + (l- a)ip{y). 
Since if is continuous it is also homogeneous and convex on M*^. □ 

4.4 Behavior of r 

Our next result says that for z G Z'', r(0, z) grows at most linearly in ||-2||oo- 

Theorem 4.2 There exist K = K{X, d) > and a > such that 

P{r{o,z)>K\\z\U < exp(-«(||^||^/^), V ^ G 
P{t{o,z) > K{\\z\\^ + n)) < exp(-an^/"'), V^GZ^nGN, 
^P{t{o,z) > K\\z\\^) < +00. 



Proof of theorem 4-2. Let K > 0, z ^ Z"^. Then write: 



P(f(o, z)> K{\\z\\oo + n)) 
< P{4K{z)>\\z\\^+n) + P{4n{o)>\\z\\^ + n) + P{A) (43) 



where 



A = {r{o,z)> K{\\z\\oo + n),4:K{z) <\\z\\oo + n,4K{o) <\\z\\oo + n} (U) 
C l^(cc,y)eB{o,(M\oo+n)/4)xB{z,(M\oo+n)/4){x ^ ?/, r{x, y) > K{\\z\\^ + n)}. 

Noting that if {x, y) G B{o, {\\z\\ao + ^)/4) x B{z, (||z||oo + n)/A) we have 

Iklloo-^ < 2||2; - ylloo < 3||2;||oo + ^ and 
^{\\z\\^ + n) = 3||2;||oo +n + 2n > 2(||x-y||oo +n), (45) 



from (44), for C2 given in Lemma 3.5 we get: 



x(^B{o,{\\z\\oo+n)/i) yeB{z,{\\z\\^+n)/A) 

{p{?>T{x,y)>2K{\\x-y\U + n),D{x,y) < (C2 + l)(||x - + n) 
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+P(x ^ y, D{x,y) > (C2 + l)(||x - y\U + n))). 



(46) 



It now follows from Lemma 3.5 part (i) that we have 



P{x ^ y, D{x, y) > (C2 + l)(||x - y\\i + n)) 

< exp{-as{\\x-y\\i + ny/'') 

< expi-asiWx -yW^+riY^'^). 



(47) 



Then, taking K large enough, by large deviation results for exponential vari- 
ables, we also have, for some > 0, 

P{3t{x, y) > 2K{\\x -y\\^+n), D{x, y) < {C2 + l){\\x - y\\, + n)) 

< P{3T{x,y) > 2K{\\x -y\\^+n),D{x,y) < {C2 + l)d{\\x - y\\^ + n)) 

< exp(-a5(||a;-?/||oo + ^)). (48) 



Hence, from (45)-(48), for some constants R and ae > we have: 



PiA) < R{\\z\\^ + n)2^exp(-«6(||2|U + n)^/^), 
which gives, by modifying the constants, 

PiA) < i?'exp(-a7(lkl|oo + ny/''). (49) 



All the statements of the Theorem now follow from (49), (43) and Lemma 
K2\ □ 

4.5 Asymptotic shape for r 

Theorem 4.3 Let e > 0, and 



At = {z eZ'^ ■.t{o,z) <t}, 
D = {x eM."^ : ifiix) < 1}. 

Then, a.s. fort large enough, 

(1 - e)tD n Z'^ c A c (1 + £)tD n z^. 



(50) 



(51) 



Remark 4.5 The set D is bounded: indeed passage times along edges are 
bounded below by passage times of exponential distributions, hence the epi- 
demic cannot propagate quicker than this first passage percolation process, 
whose passage times have exponential distribution of parameter X/{2d), and 
which, by ^ Theorem (1.15)], moves linearly following the boundary of a 
convex set. 
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In the sequel K is a fixed constant satisfying tfie conclusions of Theorem 



4.2, 7* is the Lipschitz constant of (p (see (36)) and A^^ was defined in (38) 
for any a; G Q"' \ {o}. 

Lemma 4.10 Let p > and let 5 < p/{2K). Then, for all x eQ'^\ {o}, 

y'P( sup T{kN^x,z) > kN^p) < oc, (52) 

zeBikN^x,SkN^)nZ'i 

y'P( sup T{z,kN^x) > kN^p) < oo. (53) 

z<^B{kN^x,5kN^)rM.'i 



we 



Proof of lemma \4.1(\ Let A; > 0, ;2 G B{o, 5kN.^) n Z'^. By Theorem [42 
have: 

P{r{o,z) > kN.,p) < P{t{o,z) > K\\z\\oo + [kN^p/2\) 
< exp(-a[fciV,p/2ji/'^). 

Therefore, for some constant C, 

Vp( sup r{o,z)>kN,p) < VC((5A;iV,)'^exp(-«[A;iV,p/2ji/'^) 

fc>o zeB{o,SkN.)nZ'' fc>o 

< oo. 



Now (52) follows from the translation invariance of r. The proof of (53) is 



analogous. □ 

For X = (xi, . . . , Xd) G Q'^ \ {a} and 5 > 0, we define the cone associated to 
X of amplitude 6 as 

C{x,S) = Z'^r\ (^Ut>oB{xt,St)y (54) 



Lemma 4.11 Let x e Q'^\{o}. Then for any < S' < S the set C{x, 5') \ 
Uk>oB{kNxX, 6kNx) is finite. 

The proof of this lemma is elementary and left to the reader. 



Proof of theorem 4-3 Fix e G (0, 1) and let p,6 and l be three small posi- 
tive parameters such that S < p/ {2K), whose values will be determined later. 
The set 3^ = {x G Q*^ : 1 — 2l< (f{x) < 1 — i} is a ring between two balls with 
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the same center but with a different radius, because by Lemma [4. 9[ ip is ho- 
mogeneous and positive except that ip{o) = 0. Hence the (compact) closure 
of 3^, which is recovered by balls of the same radius centered on the rational 
points of 3^, is in fact covered by a finite number of such balls. Thus there ex- 
ists a finite subset y of 3^ such that Z'' C Ux£yC{x, 6/2) (if the balls recover 
the ring, the cones associated to them recover the whole space). Hence, to 



prove the first inclusion of (51 ) it suffices to show that for any x eY and any 
sequences (tn)n>o and {zn)n>o such that t„ t oo in M"*", Zn G C{x,6/2) fl Z'^ 

< tr, a.s. for n suf- 



with 



> n and (p{zn) < (1 — £)tn, we have r(o, Zn) 



ficiently large. So, let (t„)n>o and {zn)n>o be such sequences. Using Lemma 



4.11 , let /c„ G N be such that Zn G B{knNxX, 5knNx), hence kn > Cn for some 
constant C. Since by Lemma [4.9[ ip is Lipschitz with Lipschitz constant 7*, 
write, for 7 = 7*(i: 

knNx{l - 2l) < p{knN-,x) < p{Zn) + l^knN^ < (1 " ^)^n + l^k^N^. 

Therefore 



kuNx < 



l-e 



1 - 2i - 75 



It now follows from this inequality and the subadditivity property (21) of r 
that: 



T O, z 



1 '^n ) 



< 



1 - 2i - 7(5 



f(0, knNxX) ^ u{knNxX) ^ T{knNa: X, Zj^ 



k-n-^x 



distributed, and kn > Cn), Lemmas 4.9 and 4.10 we obtain: 



Therefore, by Theorem 4.1, Lemma 4.6 (the variables u{.) are identically 



lim sup 



T O, Zr, 



tr. 



< 



l-e 



Since x eY this implies: 
lim sup 



n—^+oo 



< 



2i-75 



1-6 
-2i--ib 



p[x) + p 



I + p 



a.s. 



a.s. 



Taking t, p and 5 small enough, the right hand side is strictly less than 1 
which proves that r(o, 2;„) < tn a.s. for n sufficiently large. 



Similarly, to prove the second inclusion of (51 ) it suffices to show that for 
any x eY and any sequences t„ t oo in M"*" and Zn in C(x, 5/2) fl Z'^ such 
that p{zr/) > (1 + e)tn we have r(o, z„) > tn a.s. for n sufficiently large. As 
before, we let (t„)„>o and {zn)n>o be such sequences and we let kn E N he 
such that Zn G B{knNxX,6knNx). Then, 



knN:^{l - > piknN^x) > p{Zn) - jSknN^ > (1 + ^)tn " l^knN^ 
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Therefore, 

l + e 



knN, > -^]tn. 

\1 — L + 'JO/ 

Proceeding then as for the first inclusion, we get: 

r(0, Zn) ^ / 1 + g \ / f(0, knN^x) _ u{Zn) _ T{Zn, 

and 

hminf ^ ' > ( — 7:^('^(x) — p] a.s. 

Now, taking l, p and 5 small enough, the right hand side is strictly bigger 



than 1 and the second inclusion of (51) is proved. □ 



4.6 Asymptotic shape for the epidemic 

We can now prove our main result: 



Proof of theorem 2.1. (i) We first show that the infection grows at least 



linearly as t goes to infinity, that is, given e > 0, 

p{{Ct U it) 3 ((1 - e)tD n C°) for all t large enough) = 1. 
Since R° is finite a.s. this will follow from: 

p{{Ct U it) 3 ((1 - e)tD n (C° \ for all t large enough) = 1. (55) 



Let 2 G (1 - e)tD n {C° \ R°), then by Theorem 4.3 



t{o,z) < (1 — £/'2)t, a.s. for t large enough, (56) 

and by Lemma [43 , r(o, z) < {1 — s/2)t + u{o) + u{z). Since u{o) < oo a.s. 
we have u{o) < {e/A)t a.s. for t large enough. Hence (55) will follow if we 
show that svLp^(ztDu{z) < {e/4:)t a.s. for t large enough: To derive this, it is 
enough to show that sup^g(„_,_i)£, m(z) < (£/4)n a.s. for n(e N) large enough. 
By Remark 4.5, D is bounded, hence the number of points in (n + 1)D with 
coordinates in Z is less than C^ln + l)'^ for some constant C5. Then write 



P I sup u{z) > ^ I < C^in + 1)''P (u{o) > 
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Thus, by Lemma 4.6, J2neN Pi^^Pze{n+i)D '^i^) — ^n/A) < oo, and (55) fol- 
lows from Borel-Cantelli's Lemma. 

(a) Next we show that 

p((Ct U 6) C ((1 + e)tD n C°) for all t large enough) = 1. (57) 



If z belongs to or then r(o, z) < t hence by Lemma 4.3, r(o, z) < t for 



z G C°\ R°, which implies 2; G (1 + e)tD for t large enough by Theorem 4.3 



Since R° is finite (57) follows. 



(Hi) Finally, assuming E{\Tz\'^) < 00, we show that 

P{Ct n (1 - e)tD = for t large enough) = 1. 



(5^ 



Let z e {1 - e)tD n C°, then, by (56), t{o, z) < (1 - e/2)t if t is large 



enough. Hence, (58) follows if we show that > {e/2)T(o, z) occurs only 



for a finite number of z's. But from (57) we get that for some 5 > we 



have t{o, z) > 6\\z\\oo except for a finite number of z's. Therefore, it suffices 
to show that for any 6' > the event {T^ > 5'||2||oo} can only occur for a 
finite number of z^s. This will follow from Borel-Cantelli's Lemma once we 
prove that ^2gz<*-P(^2 — ^'Iklloo) < 00. To do so we write, since the T^'s 
are identically distributed: 

J2Pin>S'\\z\U = J2 Yl PiTz>5'n)<cY,n''-^P{To>5'n) 

n6N 2::||2||oo=" 



for some constant c, and this last series converges because Tq has a finite 
moment of order d. □ 
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